SECURITY 


3CAHON  OF  "ms  7aGE 


la.  REPORT  SECURITY  CLAS$IFlCAT( 

a 


2a  SECURITY  CLASSIFICATION  AU 


2b  DECLASSIFICATION  ,  DOWNGRA 


4  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 

Technical  Report  #4  -■ 


DOCUMEf 

i  ELECT E 


AD-A223  614 


3  OlSi.  _ .vmilABILiTY  OF  REPORT 

Approved  for  public  release. 
Distribution  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a  NAME  OF  PERFORMING  ORGANIZATION 

University  of  Minnesota 


6o.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANlZAT.ON 
(If  applicable) 

Office  of  Naval  Research 


6c  AOORESS  (City,  State,  and  ZIP  Code) 

Department  of  Chem.  Eng.  &  Mat.  Sci. 
421  Washington  -Avenue  SE 
Minneapolis,  MH  55455 


7b  AOORESS  (Ciry.  Stare,  and  ZIP  Code) 

800  N .  Quincy  Street 
Arlington,  VA  22217-5000 


9a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


8c.  AOORESS  (City.  State,  and  ZIP  Code) 

800  N.  Quincy  Street 


Bb.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

N00014-88-K-0366 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

5000 

ELEMENT  NO 

NO 

NO 

ACCESSION  NO 

i 

1 1  TITLE  (Include  Security  Classification) 

Calculation  of  Average  Koleular  Properties  During  Nonlinear,  Living  Copolymerization 


j  12  PERSONAL  AUTHOR(S) 

1  C.W.  Macosko  and  D.R.  Miller 


1 3a.  TYPE  OF  REPORT 

Interim 


16  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM  _  TO 


14  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 

1990  June  25 


LfcMtNIAKT  rwvjiAinjr^  ^  ArTuz.-fon^, 

Revised  version  of  a  paper  submitted  to  Die  Makromolekulare  Chemie,  April  1990 


COSATI  CODES'’ 


GROUP  A'  SUB-GROUP 


19.  ABSTRACT  [Continue  on  reverse  if  necessary  and  identify  by  block  number) 


Describing  network  formation  during  chainwise  polymerization  is  more  complex 
than  for  stepwise  (condensation)  systems.  The  kinetics  of  chain  growth  must  be  combined 
with  probability  laws  for  random  linking  between  chains.  This  approach  is  used  to  derive 
molecular  weight  averages,  gel  point  and  sol  fraction  during  living,  anionic 
copolymerization  of  vinyl  with  divinyi  monomers.  Comparison  of  predicted  molecular 
weight  and  gel  point  to  experimental  values  from  Worsfold  (1970),  Lutz,  Beinert  and 
Rempp  (1982)  and  others  indicate  that  a  significant  amount  of  intramolecular  linking  or 
cyclization  occurs  in  this  system*. 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

IS  UNCLASSIFIEDAJNUMITED  □  SAME  AS  RRT  □  DT1C  USERS  I  Unclassified  _ 


22*  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b.  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 

Dr.  Joanne  Milliken/Dr.  Kenneth  Wynne  202-696-4409 


00  FORM  1473, 84  MAR  83  APRM/tion  maybe  u»«d  until  «*n*u«ed.  security  n  aerification  of  this  page 

All  other  editions  are  obsoiRt*. 


Revised  4/13/90 


CALCULATION  OF  AVERAGE  MOLECULAR  PROPERTIES  DURING  NONLINEAR, 

LIVING  COPOLYMERIZATION 


Dedicated  to  Paul  Rempp  in  honor  of  his  60th  birthday 


by 

Christopher  W.  Macosko 

Department  of  Chemical  Engineering  and  Materials  Science 
University  of  Minnesota 
Minneapolis,  MN  55455 

and 


Douglas  R.  Miller 

Department  of  Operations  Research  and  Applied  Statistics 
George  Mason  University 
Fairfax,  VA  22030 


Summary 

Describing  network  formation  during  chainwise  polymerization  is  more  complex  than  for 
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predicted  molecular  weight  and  gel  point  to  experimental  values  from  Worsfold  (1970),  Lutz,  Beinert  and 
Rempp  (1982)  and  others  indicate  that  a  significant  amount  of  intramolecular  linking  or  cyclization  occurs 
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Submitted  to  Die  Makromolekulare  Chemie,  Revised  March  1990. 


2 


INTRODUCTION 

The  development  of  molecular  structure  during  nonlinear  polymerization  can  be  modelled  with 
reasonable  accuracy  for  stepwise  (condensation)  growth.  For  the  ideal  case  all  functional  groups 
combine  randomly,  independent  of  the  state  of  neighboring  groups  or  size  of  the  molecule  to  which  they 
are  attached.  Intramolecular  reactions  do  not  occur  between  finite  species.  With  these  assumptions 
nonlinear  stepwise  polymerization  can  be  modelled  by  branching  theory.  The  branching  probability 
model  can  be  solved  by  combinatorial  methods1,  probability  generating  functions2  or  directly3  to  calculate 
parameters  such  as  molecular  weight  averages,  gel  point,  soluble  fraction  and  crosslink  density.1'6  Time 
only  enters  the  problem  through  a  single  equation  relating  the  conversion  of  functional  groups  to  the 
reaction  rate. 

A  number  of  real,  bulk  condensation  systems  do  polymerize  nearly  ideally.  Experimental 
measures  of  weight  and  number  average,  gel  point  and  sol  fraction  agree  well  with  branching  model 
results  (e.g.  ref.  7-10).  When  functional  groups  are  not  equally  reactive,  kinetic  equations  to  describe  the 
various  linking  probabilities  for  some  cases  can  be  solved  analytically  and  incorporated  into  the  branching 
model.4*11-12  Intramolecular  cyclization  generally  is  only  important  in  dilute  systems.  Even  for  this  case 
introduction  of  kinetic  equations  for  the  smallest  size  cycles  provides  useful  results.13*14 

Nonlinear  polymerization  by  chainwise  addition  is  more  difficult  to  model.  In  general  chains  do 

not  grow  randomly  so  kinetic  equations  must  be  solved  for  the  length  distribution  of  the  primary  chains  in 

the  network  and  then  a  probability  model  devised  to  combine  these  chains  randomly  into  a  network.  For 

the  special  case  of  slow  chain  initiation  compared  to  propagation  (i.e.  free  radical  polymerization)  with 

chain  termination  by  disproportionation  and  equal  reactivity  of  the  different  vinyl  groups  Flory1  showed 

that  the  random  branching  model  can  be  applied.  However,  when  there  is  termination  by  coupling  the 

direction  of  propagation  along  each  chain  must  be  considered.15  Unequal  reactivity  can  be  treated  with 

kinetic  equations  for  the  various  linking  probabilities  in  a  manner  similar  to  that  used  for  condensation 

systems.  Cyclization  is  very  important  in  nonlinear  free  radical  polymerization  (e.g.  ref  16).  Long 

chains  with  many  pendent  functional  groups  form  early  in  the  reaction.  Thus  the  network  is  typically 

_ _  deal 

formed  in  dilute  solution.  /  o,lc\  and/or 

|  co»y  )  Dlst  j  -Peoiai 

iUtaa _  I 
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In  this  paper  we  are  interested  in  "living"  anionic  chain  copolymerization.  Linear  anionic 
copolymerization  has  been  treated  by  solving  the  kinetic  equations17  and  by  Monte  Carlo  simulation.18 
The  gel  point  for  nonlinear  copolymerization  with  equal  reactivity  was  derived  by  Eschwey  and 
Bure  hard.1 9  Bokare  and  Gandhi20  treated  combined  stepwise  and  living  polymerization.  Dusekand 
Somvarsky21  modelled  molecular  weight  up  to  the  gel  point  for  non-linear  anionic  homopolymerization, 
but  it  appears  that  no  one  has  modelled  the  general  problem  of  nonlinear  copolymerization  with  unequal 
reactivity.  Nor  has  there  been  much  comparison  to  experimental  data. 

The  anionic  homopolymerization  of  diepoxides  and  diisocyanates  have  been  used  to  make 
network  polymers.20*22  There  has  also  been  considerable  interest  in  making  star  molecules23  and 
controlled  microgel  particles24  via  anionic  polymerization  of  divinyl  monomers.19*23-26  These  studies 
have  examined  what  reactivity  ratios  and  what  vinyl  to  divinyl  ratios  will  prevent  gelation  in  anionic 
polymerization  of  divinyl  monomers:  i.e.  how  can  a  controlled  amount  of  branching  be  introduced? 
Several  studies  indicate  that  cyclization  is  important  in  these  systems.19*29 

In  this  paper  we  treat  "living"  anionic  copolymerization  using  the  kinetic-Markovian  approach 
described  recently.27  A  kinetic  model  describes  the  growth  of  the  primary  chains  which  are  connected 
randomly  via  a  branching  model  to  give  a  network  polymer.  The  model  includes  the  possibility  of 
unequal  rates  of  addition  of  the  different  types  of  vinyl  groups  and  we  indicate  how  other  complications 
such  as  slow  initiation  can  be  treated.  We  derive  relations  for  number  and  weight  average  molecular 
weight  up  to  the  gel  point  and  weight  fraction  solubles  beyond  the  gel  point.  These  relations  are 
compared  to  experimental  results  from  the  literature. 

We  consider  a  system  consisting  of  two  types  of  reactants,  A2  and  A4,  and  an  initiator  I.  We 
use  the  notation  A2  to  denote  a  monomer  with  one  reactive  site  which  when  activated  acts  as  a  2- 
functional  unit  In  this  paper  we  compare  data  for  anionic  polymerization  of  vinyl  monomers.  Therefore 
the  reactive  site  is  C  =  C  but  it  could  also  be  an  isocyanate,  epoxide  or  other  ring.  The  notation  A4 
denotes  a  monomer  with  two  reactive  sites  so  that  when  both  are  activated  the  monomer  acts  as  a  4- 
functional  unit  Divinyl  benzene  is  an  example  of  an  A4  monomer.  The  different  units  have  masses 
denoted: 
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Ma2»  Ma4>  and  Mi 
Their  initial  concentrations  are 
[A2]0,  [A4]o,  and[I]o 
respectively. 

THE  IDEAL  CASE 

Before  we  construct  the  general  model  it  is  useful  to  consider  the  simplified,  ideal  case.  If  we 

assume 

1 .  equal  reactivity  between  all  reactants:  initiator,  A2  and  both  vinyls  on  the  A4 

2.  no  substitution  effect,  all  groups  react  independently  throughout  the  polymerization  and 

3 .  no  cycles  form  in  finite  species. 

These  are  the  same  assumptions  as  for  the  ideal  stepwise  case  mentioned  above.  In  addition,  we  will 
make  three  additional  assumptions: 

4 .  The  rate  of  activation  of  the  initiator  is  very  fast  compared  to  the  growth  rate  and  is 
assumed  to  be  instantaneous. 

5 .  The  initiation  rate  is  the  same  as  the  growth  rate. 

6.  The  mass  of  the  initiator  is  negligible  and  is  assumed  to  be  zero. 

Recently  Gupta  and  Macosko30  have  solved  the  case  of  Af  polymerization  with  initiators  Ig.  The  general 
case  of  any  ratio  of  initiation  to  growth  has  been  solved. 

With  these  assumptions  we  can  separate  the  living,  anionic  chain  growth  from  crosslinking  and 
treat  crosslinking  as  a  completely  random  process.  Initially  we  will  ignore  the  connection  between  the 
two  vinyl  groups  on  the  A4S  and  solve  the  kinetic  equations  for  the  linear  chain  length  distribution. 

Since  all  species  are  equally  reactive  and  have  the  same  mass,  we  do  not  need  to  distinguish  between  A2 
and  A4  until  the  crosslinking  step.  Let  [A2]0+  2[A4]0  =  [MJo,  the  initial  concentration  of  monomers.  The 

initiation  and  growth  reactions  can  be  written  as 


k 

I  +  M  -»Pi 


(1) 
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Pn  +  M  -»  Pn  +  i  ,  n  =  1(2,3, ... 

where  Pn  is  a  living,  linear.chain,  n  monomer  units  long.  The  rate  equations  become 


^U-k[M][H 


^ll  =  k[M][I]-k[M][P1] 


®^  =  k[M][Pn.i]-k[M][Pn] 


for  n  =  2,3,. 


Time  and  the  rate  constant,  k,  can  be  eliminated  from  these  equations  by  considering  the  overall  decrease 


of  monomer 


^U-ktMinio 


which  when  divided  into  the  above  gives 


dm  _  [i] 

d[M]  [Ho 


dPL_JEL+mil 
dM-[Ho  [I]0 


_  ipqji  +  jai 

d[M]  ~  [I]0  +  [I]0 


for  n  =  2,3,. 


whose  solution  is 


m  =  m0e'4 

[Pi]  =  chohc'^ 

[Pn]  =  [I]#V>Vn!  for  n  =  2,3,... 
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where 

[M]0  -  [M]t 
“ '  (Ho 


or  in  terms  of  the  conversion,  a: 


reacted  vinyl  groups 
initiator 


M-  =  a[My[I]0  (8) 

This  number-distribution  (Equations  6  and  7)  is  the  well-known  Poisson  distribution,  just  the 
result  that  we  expect  for  anionic  living  polymerization.28  The  mean  of  the  distribution  is  |i,  the  number 
average  chain  length.  The  actual  chains  may  consists  of  0,  1,2, ...  monomer  units.  (A  chain  of  0  length 
corresponds  to  an  unreacted,  activated  initiator.)  Furthermore,  because  A2's  and  A^s  join  living  chains 
independently  and  at  equal  rates,  the  fraction  of  crosslinkable  monomers  (A4S)  in  the  chain  is 

P  =  2([A4]0/([A2]0  +  2[A4]o)  =  2[A4]o/[M]0  (9) 

and  the  number  of  A4S  in  a  random  chain  has  a  Poisson  number-distribution  with  mean  pp.  Similarly 
the  number  of  A2S  in  a  chain  has  a  Poisson  number-distribution  with  mean  (l-p)p,  independent  of  the 
number  of  A4S. 

Now  let  us  put  back  in  the  links  (between  the  vinyls)  on  the  A4  monomers.  These  will  become 
crosslinks  between  chains.  Under  our  idealizing  assumptions  it  is  possible  to  derive  formulae  for  Mw 
and  etc  at  gel  using  the  recursive  method.3,6  We  first  randomly  pick  a  root  by  mass;  it  will  be  an  A4  or 
A2  with  probabilities 

_ [A4]0  (1  -  (1  -  oQ2)  MA4 

WA4  [A4]0  (1  -  (1  -  a)2)  MA4  +  [A2]q  a  MA2 

_ [A2q]  Ct  MA2 _ 

WA2*  [A4]0(l  -  (1  -  a)2)  MA4+  [A2]0a  MA2 


(10) 
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respectively.  Note  that  [A2]Q  a  is  the  concentration  of  reacted  A2's  and  [A4]0  (1  -  (1  -  a)2)  is  the 

concentration  of  A4*s  with  at  least  one  reacted  vinyl  group.  The  denominator  above  is  the  usual  polymer 

"yield,”  the  mass  after  unreacted  monomers  have  been  removed. 

Let  the  additional  average  weight  of  the  branched  chain  through  a  reacted  A-unit  be  denoted  as 
E(WA ,).  Using  the  recursive  method  and  the  fact  that  the  average  number  of  additional  A2's  and  A^s  on 

a  chain  are  (l-p)p.  and  pp,  respectively,  we  have 


E(WA2)  =  MA2  +  E(W* ,) 


(12) 


If  the  chosen  root  is  an  A4,  we  know  that  there  is  at  least  one  chain  through  the  unit  and  possibly  two. 

For  a  randomly  chosen  A4  we  get  0,  1,  or  2  chains  with  probabilities  (1  -  a)2,  2a(l  -  a),  or  a2, 
respectively;  thus  conditioning  on  the  A4  being  part  of  the  yield  we  have  1  or  2  chains  with  probabilities 
2a(l  -  a)/(l  -  (1  -  a)2  or  a2/(l  -  (1  -  a)2,  respectively,  so  the  expected  number  of  chains  equals 
[2a(l  -  a)  +  2a2]/([l  -  (1  -  a)2]  =  2a/(2a  -  a2)  =  2/(2  -  a).  Thus 

EOV^)  =  +  [2/(2  -  a)]  E(WA  t)  ( 1 3) 

Combining  equation  10  through  13  gives 

Eiw  =  wA2  E(WA2 )  +  w^  EfWyj^)  14) 

Note  that  this  is  the  weight  average  molecular  weight  of  the  yield. 
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In  order  to  compute  the  weight  average  molecular  weight  of  the  entire  species  the  above  analysis 
must  be  modified,  using  the  following  equations 


WA4  = 


wA2  = 


[A4]q  Mm 


[a4]g  mA4  +  [a2]0  mA2 


(10') 


[A2]n  MA2 


[a4]0  mA4  +  [A2]q  mA2 


E(WA2)  =  MA2  +  oE(W^) 


(12’) 


E(WA4)  =  MA4  +  2aE(W^ut)  (13’) 

The  resulting  weight  average  molecular  weight  using  equations  10’,  11, 12’,  13’,  and  14  will  be  for  the 
entire  species  rather  than  for  the  yield. 

From  equation  1 1,  we  see  that  the  average  weight  diverges  and  we  have  gelation  at  conversion 
etc  =  1/pM-  (15) 

Combining  equations  8, 9  and  15  gives 


oc  =  VtI]o/2[A4]0 


(16) 


Note  that  2[A4]y[I]o  equals  the  ratio  of  the  number  of  potentially  crosslinkable  chain  segments  to  the 

number  of  chains;  this  number  must  be  greater  than  1  in  order  for  gelation  to  occur. 

Now  we  restrict  our  attendon  to  systems  consisting  of  only  A4  units.  The  weight  average 
molecular  weight  of  such  an  anionic  polymerization  from  equations  10  through  14  becomes 
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Bw  =  E(Wa4)  -  MA4  +  MA4  *  (77^  )T 


^MA4 


a  (i 


Mw  =  M 


+ 


- 2ji - 1 

(2  -  a)  (1  -  ap)J 


(17) 


This  is  the  weight  average  molecular  weight  of  the  yield.  If  we  consider  the  entire  species, 
including  unreacted  Aa's  we  get 


Mw=E(WA4)=  Mm  +  2aE(W*t) 


=  +  2  a 


MM 

1-aji. 


=  M 


M 


fi+-22M 

V  l*«nj 


=  M 


A4 


rl  +  ap.' 

J  -  an 


(18) 


the  weight  average  molecular  weight  of  the  entire  species. 


Comparison  to  Step  and  Chain  Addition 

It  is  of  interest  to  compare  these  weight  averages  (equations  17  and  18)  to  the  averages  for  free- 
radical  chain  addition  homopolymerization  and  stepwise  homopolymerization,  respectively. 

For  stepwise  homopolymerization1'3 


AfU  +a) 

1  -  a  (f  -  1) 


(19) 


which  simplifies  to 
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MA4(1  +  a) 

1  -  3a 


for  tetrafunctional  monomers 


and  the  gel  point  is 


Ideal  free  radical  homopolymerization  of  multivinyl  monomers  leads  to  (eqs.  77, 79  pg.  204, 
ref.  3  with  MA2  =  0  and  aj-  =  1) 


Mw  = 


MAf(l+q) 

1-q  [1  +  (f-  2)a] 


where  q  is  the  probability  that  a  growing  chain  will  add  one  more  monomer  unit  rather  than  terminate  by 
disproportionation.  The  gel  conversion  is 


1-q 


Figure  1  compares  Mw  for  the  three  polymerization  mechanisms.  For  typical  conditions  we  see 
that  living  polymerization  shows  characteristics  of  both  the  step  and  chain  mechanisms  for  molecular 
weight  growth.  Like  the  stepwise  case,  the  initial  molecular  weight  of  the  living  chains  is  that  of  the 
monomer  (since  we  ignore  initiator  weight).  But  as  the  A4  monomers  add  to  the  living  end  the 
functionality  of  these  chains  increases  rapidly  and  they  begin  to  branch.  At  the  gel  point  the  weight 
average  functionality  of  each  living  chain  is  10. 

The  ideal  chainwise  mechanism  produces  long  chains  of  A4  in  the  first  instant  of 
polymerization.  Their  high  weight  average  length,  L*  =  100,  means  that  they  have,  on  average,  100 
functional  groups  per  chain.  It  takes  very  little  further  reaction  for  these  groups  to  connect  the  few  chains 
around  into  a  network. 
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For  all  three  mechanisms  the  gel  point  is  determined  by  the  weight  average  functionality  at  that 
point.  This  can  be  shown  by  comparing  the  three  gel  point  equations. 


living 

Zi 

n 

eq. 

(15) 

step 

JP 

n 

t— * 

i 

*— » 

eq. 

(20) 

chain 

etc  =  l/[(f— 2)  q/(l-q)J 

eq. 

(22) 

If  we  consider  the  effective  functionality  of  all  three  systems,  we  see  that  it  is  exacdy  the  weight 
average  chain  length,  L  w  (because  every  segment  is  a  potential  crosslink).  The  weight  average  length  for 

the  living  anionic  is  (from  eq.  7) 

rw=ii+i  (23) 

The  result,  with  eq.  15,  gives 

Oc=l/(Lw-l)  (24) 

which  is  the  same  form  as  for  the  stepwise  case. 

In  the  chain  addition  case,  assuming  a  most-probable  distribution  ,  we  have6  L  n  =  1/(1  -q) 
and  L  w  =  2  L  n  -  1.  This  gives 

Lw  =  (l+q)/(l-q)  (25) 

which  when  substituted  into  eq.  22  with  f  =  4  gives  eq.  24.  So  all  the  gel  point  equations  are  of  the  same 
form.  This  is  illustrated  in  Figure  1  where  L  w  is  set  at  4  for  the  living  (2[A4]q/[I]o  =  9),  and  for  the 

chain  systems  (q  =  0.6).  Although  gel  points  are  matched  we  still  see  differences  in  the  shape  of  the 
molecular  weight  rise  due  to  each  mechanism. 
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Instantaneous  Initiation 

In  deriving  the  Poisson  distribution  above  we  assumed  initiation  and  growth  reactions 
proceeded  at  the  same  rate.  However,  this  is  not  generally  the  case.  Typically  initiation  is  faster  than 

propagation  and  often  chains  are  preinitiated.  A  small  amount  of  monomer  is  added  to  the  initiator 
([MJj  -  2[I]o)  and  then  after  the  chains  have  all  started  growing  the  rest  of  the  monomer  is  put  in.  Thus  it 

is  important  to  determine  the  role  of  initiation  rate.  Let  us  assume  initiation  is  instantaneous.  Then  at  time 
=  0*,  [Ho  =  0  and  [Pi]  =  [I]o  •  Applying  these  initial  conditions  to  the  above  derivation  gives  a  slightly 
different  distribution 


Pn  =  [I]0v"-ie-v/(n- 1)!  n=  1,2,3,... 

(26) 

where 

[M]0+  •  [M]t  [M]0-[I]0-[M]t 
v  ■  [1)0  ■  [Do 

•H-  l 

This  can  be  called  a  shifted  or  alternate  Poisson  distribution.28 

The  weight  average  chain  length  for  this  distribution  is  28 

(27) 

Lw  =  (l  +  3v  +  v2)/(1  +  v) 

(28) 

L  w  =  M-  +  1  •  1/jJ- 

(29) 

This  result  can  be  substituted  into  equation  25. 

The  gel  point  becomes 

Oc  =  { [Il0/2[A4]0  +  ([I]q/2[A4]0)2  }  1/2 

(30) 
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This  result  converges  to  eq.  16  for  long  chains.  This  is  illustrated  in  Table  1.  We  see  from  the  table  that 
for  chains  with  10  or  more  divinyl  units  incorporated  at  complete  conversion  there  is  essentially  no 
difference  between  equal  and  fast  initiation.  Dusek  and  Somvarsky21  have  solved  homopolymerization 
with  arbitrary  initiation  rate.  Some  of  their  results  are  included  in  Table  1.  We  see  that  when  the 
initiation  rate  becomes  10  times  the  propagation  rate  it  is  essentially  instantaneous.  We  will  assume  fast 
initiation  in  the  remainder  of  this  paper. 

UNEQUAL  REACTIVITY 

In  anionic  copolymerization  the  monomers  are  generally  not  equally  reactive.  The  A2,  the  first 
vinyl  group  on  an  A4  and  pendent  vinyl  groups  add  to  the  growing  chain  at  unequal  rates.  To  model  this 
problem  we  make  the  following  assumptions: 

1 .  initiation  reaction  is  very  fast  compared  to  the  chain-growth  reaction,  so  that  the  initiation  can 
be  assumed  to  occur  instantaneously  at  time  t  =  0. 

2 .  second  order  kinetics  during  chain  growth. 

3 .  rate  of  a  unit  joining  one  of  the  living  chains  depends  solely  on  the  type  of  unit  and  not  the 
type  of  unit  on  the  living  end.  This  assumption  appears  valid  for  a  number  of  anionic 

copolymerizations.22.29 

4 .  units  join  chains  independently  of  one  another. 

5 .  no  intramolecular  reactions,  i.e.  no  cycles  in  the  finite  species. 

We  shall  explicitly  treat  the  case  where  the  initiator  creates  a  chain  with  one  live  end.  The  case  where  a 
chain  has  two  live  ends  is  virtually  identical;  the  chains  grow  twice  as  fast  because  there  are  double  the 
number  of  chain  ends.  Each  chain  also  has  two  initiated  units  instead  of  just  one. 

We  take  a  hybrid  modelling  approach  similar  to  those  described  by  Miller  and  Macosko28.  The 
kinetic  aspects  are  modelled  separately,  yielding  superspecies;  these  superspecies  are  then  combined 
randomly  into  molecules  and  networks  using  a  probability  model.  This  system  is  an  especially  interesting 
modelling  exercise  because  it  breaks  into  four  conceptual  stages:  the  initiation  stage,  the  chain  growth 
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kinetic  stage,  the  chain  growth  random  combination  stage,  and  the  chain  crosslinking  stage.  At  the  end  of 
the  paper  we  show  some  calculations  and  compare  to  some  experimental  results  in  the  literature. 

Initiation  Stage 

During  the  initiation  stage  the  initiator  activates  sites  on  the  A2  and  A4  monomers.  The 
unreacted  monomers  can  be  denoted  as 


Here  styrene  and  m-divinyl  benzene  have  been  used  to  illustrate  our  notation.  A  short¬ 
hand  notation  will  also  prove  useful  for  illustrating  the  growing  chain  structure 


o 


o 


These  unreacted  monomers  greatly  outnumber  the  initiator  so  after  the  initiation  stage  most  of  them  will 
remain  but  some  of  them  will  be  activated  into  the  following  species 


where  denotes  the  live  end  of  a  chain. 
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We  can  use  the  following  short-hand  notation  for  these  species 


i  i 


We  assume  that  the  initiator  attacks  all  of  the  double  bond  sites  with  equal  probability  and 
independently.  We  also  assume  a  virtually  instantaneous  initiation  stage,  so  that  initiation  is  complete  by 
time  t=0+.  This  can  be  forced  by  preinitiating  a  small  amount  of  monomer  slightly  in  excess  of  [I]  then 
adding  the  remainder. 

When  initiation  is  complete,  the  proportion  of  initiated  double  bonds  equals 

... _ 00 _ 

q>  [A2]q  +  2[A4]0 

Since  the  initiator  attacks  all  double  bonds  with  equal  probability  and  independendy,  the  concentration  of 

the  various  species  t  =  0+,  are 

[A20+  =  (l-qi)[A2l0 

[AiV  =  <litA2Jo 

[AjV  =(l-qi)2[A4]0  (32) 

[AjfV  =  2qi(l-qi)[A4]0 
[aJv  =  <if  [A4J0 

It  is  possible  that  the  initiator  attacks  sites  (i.e.  double  bonds)  with  different  rates  depending  on  whether 
they  are  on  A2's  or  A4S  and  whether  the  other  site  on  an  A4  has  already  been  initiated.  For  such  a  case 
the  concentrations  of  initiated  species  will  differ  from  the  above;  they  can  be  computed  if  a  complete 
description  of  the  kinetics  is  available. 
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Chain-Growth-Kinetics  Stage 

During  the  chain  growth  stage  unreacted  sites  seek  out  live  ends  of  chains  and  attach  themselves 
to  the  chain.  A  given  site  is  equally  likely  to  join  any  of  the  live  chains,  subject  to  the  restriction  of  no 
intra-molecular  reactions.  Different  vinyl  groups  (on  A2,  on  unreacted  A4,  and  on  A4  with  one 
previously  activated  double  bond)  react  at  different  kinetic  rates:  rlt  r2,  and  r3  respectively.  We  assume 
second  order  kinetics.  Since  we  assume  fast  initiation  and  no  termination,  the  number  of  living  chains  is 
constant,  [I]0. 

During  the  growth  stage  of  the  reaction,  previously  unactivated  double  bonds  (which  we  have 
denoted  as  "o")  are  broken  and  the  site  attaches  to  the  living  chain  and  becomes  the  new  live  end.  We 
shall  denote  this  site  with  a  "g",  to  signify  the  "growth"  stage.  During  the  growth  stage  we  get  some  new 
species,  which  we  shall  denote  as  follows: 


g 

s 

a5°=  — 

c 

4  - 

o 

i 

g 

□ 

L  Ag«=  \ 

I 

g 

4  - 

g 

The  kinetic  equations  governing  the  concentrations  of  these  and  the  initial  species  over  time  are: 

a  [a?j,  - 

alAf),  =  +  r,[I]0[A°J, 

5  [A“],  *  -  2r2[I]0IA‘J>),  (33) 


5  [A?],  >  +2r2[n0[A“lt  -rj[I]0[A?]t 
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$[*?].  -  +  r3[I]0[A*o]t 
£[*?].  *  ^jIIloIA®], 

$  [A*],  =  +  rjIDotA”], 

These  kinetic  equations  can  be  solved  by  integration.  The  solution,  for  t  >  0,  is: 

[Ajt  =  [Aj0+  exp  (-r,[I]0t) 

*  [Ajt 

[A‘]t  =  [A;]0+ 

[A~],  =  [AjV  exp  (-2r2[I)0t)  (34) 

r a ^  J  [A~]o+  7rr2rT [cxp (-^rao*) - cxp c^moO] 

[A4  lt  =  i  J  1 

[or  [A°4°]0+  2r2[IJ0  exp  (-r3[I]0t),  if  r3  -  2r2  =  0 
(A“]t=[A7)0+  -[A~]t-[A*°)t 
[A4°]t  =  [AjV  exp  (-r3[I]0t) 

[a‘‘),  =  [A”V  -  [Aj0], 

[a“],  -  [A4V 

Returning  to  our  assumption  of  fast  initiation,  if  we  define  conversion  to  be  the  proportion  of 
activated  double  bonds  at  time  t,  we  get 

[A^V  [A‘]t+  [A*;],  +  [A4°]t  +  2[a“),  +  2(A'4‘]t  +  2[A“], 
m  *  (A2]0  *  2[A4J0 


(35) 
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If  we  define  yield  as  the  percentage  mass  of  the  reacted  species,  we  get 


Yield  = 


W 


Ma2[A2]0  +  MA4[A4]0  +  Mi[I]0 


(36) 


where  W  is  the  total  mass  of  polymer  (We  no  longer  neglect  the  mass  of  the  initiator.  Mi.): 


W  =  MA2  ([A*]t  +  [A  jo  +  MA4([Ag4\  +  [A4]t  +  [A“]t  +  [A*]t  +  [A«])+  Mr[I]c 


(37) 


The  number  average  molecular  weight  can  be  calculated  from  its  definition,  the  weight  per 
polymer  molecule: 


Mjj  =  W/[N] 


(38) 


The  total  number  of  chains  is  [I]0  less  the  number  of  crosslinks: 


[N]  =  [Iq]  -  [X] 


(39) 


[X]  is  just  the  concentration  of  doubly  reacted  A4  monomers: 


[X]  =  [a“]  +  [a;V(a“] 


(39) 


Thus  eq.  38  becomes 


Mn  = 


W 


[I]0  -  [A**]  -  [A1/!  -  [a;1] 


(38) 


19 


Chain>Growth>Random-Combination  Stage 

At  this  point  we  have  all  the  concentrations  of  species  from  our  kinetic  analysis.  We  now  make 
up  our  linear  chains  by  randomly  combining  our  two  monomers.  In  the  next  section  we  allow  these 
chains  to  crosslink  by  making  the  connections  between  them  through  the  reacted  A4  groups. 

Each  chain  consists  of  one  initiated  segment,  denoted  by  "i",  and  a  random  number  of 
segments  which  attached  themselves  during  the  growth  phase,  denoted  by  "g": 


8 


The  initial  segment,  "i",  may  belong  to  one  of  four  different  species 


The  growth  segments,  "g's",  may  be  pan  of 

a8  a  8°  A  *8  A  88 

/V,  ,  ,  /\4  ,  /\4 

At  time  t,  the  mean  number  of  each  "g"  segment  per  chain  are,  respectively, 
nl  -  [A|i,/ra0 

»**0  -  lAjVm,, 

H,  =  [A‘V  [Ilo  (40) 

H?  -  2(A“]t/[n0 

where  [I]o  equals  the  concentration  of  chains.  Thus  the  number  average  length  of  the  chains  is 
En(L)  =  l+p28  +  it*%p‘g  +  n84* 


(41) 
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For  any  randomly  (by  number)  chosen  chain  the  distribution  of  the  number  of  a  particular  type  of  "g" 
segment  is  Poisson.  For  example,  if  L|  equals  the  number  of  A|  units  in  the  chain,  for  £  =  0, 

1,2 . 


Pn(Lj  =  £)  =exp(-pJ)0i*)*/£! 


(42) 


The  Poisson  distribution  follows  from  the  well-known  Poisson  approximation  to  the  Binomial 

distribution:  suppose  we  pick  a  chain  randomly  (by  number)  by  randomly  choosing  an  "i"-segmenL 
There  are  [A|]t  g-segments  of  type  A2  in  one  unit  of  volume.  In  this  unit  of  volume  there  are  [I]0  chains 

to  which  the  A|  may  belong,  all  being  equally  likely.  Thus  a  given  "g"-segment  chooses  our  chain  with 

probability  l/[I]o-  This  random  choice  is  repeated  independently  by  each  of  the  [A|]t  segments.  The  total 

number  that  join  our  chain  will  have  a  Binomial  distribution  with  number  of  trials  equal  to  [A|]t  and 

probability  of  success  equal  to  l/[I]o,  yielding  a  mean  of  [Af]t  /[%  It  is  well  known  that,  when  the 


number  of  trials  is  large  and  the  probability  of  success  is  small,  the  Poisson  distribution  is  an  excellent 


approximation2^.  So,  the  length  of  a  random  chain  is 


L  =  1  +  L| +  L*u  +  L“‘  + 


(43) 


where  the  number  distribution  corresponds  to  the  sum  of  independent  Poisson  random  variables  with 
means  pf,  p8°,  p^g  and  pj8,  respectively. 

If  we  allocate  one  half  the  mass  of  A^,  A^8,  and  A88  units  to  each  chain  to  which  they  belong 

then  we  can  compute  the  average  chain  weights.  Let 
LA  s  initiated  A  unit 


P  ( IA  =  A’ )  =  [ A12]t  /  [I]o 
P(IA  =  a‘°)  =  [A;°]t/[I]o 
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P  (IA  =  A j)  =  [A'4*]|  /  [I]0  (44) 

P(IA-Af)  =  2[A’‘]t/[I]0 


if  Mia  =  mass  of  initiator  segment  on  chain,  then 


EnCM^)  =  Mj  +  P(IA  =  a‘)MA2  +  P(IA  =  a‘°)  MA4  +  P(IA  =  A*  or  a“  )  /  2  (45) 

Let 

GS  =  totality  of  growth  segments  on  chain 
and 

Mcs  =  mass  of  all  growth  segments  on  chain 

then 

En(Mos)  =  ^  MA2  +  <  +  MA4  +  (tiif  +  tig48)  MA4  /  2  (46) 

and  since  the  chain  consists  of  initiator  and  growth  segments 


E„(CW)=En(MiA)+En(MGs) 


(47) 


is  the  number-average  chain  weight,  which  should  also  equal  "Yield”  divided  by  [I]0. 

Now  consider  the  weight  average  chain  weight:  a  unit  of  mass  is  chosen  randomly  and 
designated  as  a  root,  R.  The  probabilities  of  getting  various  species  as  a  root  are 

P(R  =  A°)  =  M^A^/TM 
P(R  =  A2)  =  ( +  Mj )  [Aj]t / TM 
P  (  R  =  A8 )  =  M^A^/TM 
P(R  =  A°°)  =  MA4[A°°],/TM 
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P  (R  =  A“ )  =  (Ma4  +  Mj )  [A‘°]t/TM  (48) 

P(R  =  a;°)  =  M^A^/TM 

P(R  =  Aj)  =  (  +  Mj  )  (A*^]t  /  TM 

P(R  =  A“)  =  (MA<  +  2MI)[A“]t/TM 

P(R  =  Af)  =  MA4[Af]t/TM 

where 

TM  =Ma2[A2]0  +  Ma4[A4]0  +  Mi[I]o  (49) 

is  the  concentration  of  total  mass.  If  the  randomly  chosen  root  is  an  A^  or  an  A°°,  we  don't  get  a  chain. 

Because  of  this  we  renormalize  these  probabilities  to  just  include  the  yield  of  the  reaction: 

P'(R  =  a!j)  =  P(R  =  A‘2)/ Yield 

P'(R  =  A®)  =  P  ( R  =  A8)  /  Yield 

P'  (R  =  A4°)  =  P  (  R  =  A^°)  /  Yield 

P'(R  =  A8°)  =  P(R  =  Ag4°)/ Yield  (50) 

P'  (  R  =  A4)  =  P(R  =  A48)/ Yield 
P'  (  R  =  A“)  =  P(R  =  a“) /Yield 
P'  (R  =  Af)  =  P (  R  =  Af ) /  Yield 

For  each  root  chosen  we  get  the  root  plus  either  1  or  2  chains;  each  chain  will  consist  of  a  random 
additional  number  of  "g"-segments,  and  an  "i"-segment  if  the  root  does  not  include  the  initial  segment. 
Let 

XGS  =  Extra  growth  segments  of  chain 
XIS  =  Extra  initiator  segment  of  chain 
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MXGS  =  mass  of  XGS 
MXIS  =  mass  of  XIS 
MRC  =  mass  of  root  part  of  chain 
it  follows  that  the  total  mass  of  the  chain  is 

CW  =  MRC  +  MXIS  +  MXGS  (51) 

and  the  weight  average  mass  is 

E^CW)  =  Ew  (MRC)  +  Ew  (MXIS)  +  Ew  (MXGS)  (52) 

These  terms  can  be  evaluated  as  follows  by  conditioning  over  the  seven  different  roots: 

Ew  (MRC)  =  Sew  (MRC  I  R)  P*  (R) 

=  P'  (R  =  Aj)  (  M^  +  Mj )  +  P'  (R  =  A* )  M^ 

+  F(R-aJ,)(MA4  +  Mi)  +  F(R«aJ9)MA4  (53) 

+  P,(R  =  Aj)(M^  +  MI)/2  +  F(R*Af)M^/2 
+  P*(R  =  a“)(Ma<  +  2Mi)/2. 


Ew  (MXIS)  =£  Ew  (MXIS  I  R)  P’  (R) 

=  P'(R  =  A5)E(MIA)  +  P(R  =  Ag4°)E(MIA)  (54) 

+  P(R  =  A*)E(MIA)/2  +  F(R  =  Af)E(MIA) 

The  extra  growth  portion  of  the  chain  is  the  same  regardless  of  whether  the  chain  is  chosen  by  number  or 
weight,  thus 

Ew  (MXGS)  =  En(MGS)  (55) 

completing  the  formula  for  Ew  (CW). 
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Crosslinking-of-Chains  Stage 

Finally  we  consider  the  chains  of  the  previous  section  crosslinked  to  form  large  molecules  and 
possibly  even  a  network.  The  "crosslink"  units  are  the  species  A4*‘,  A4'8,  and  A4ss. 

It  is  tempting  to  aggregate  the  three  types  of  crosslinks  into  one  general  crosslink  thereby 
forcing  the  system  into  the  framework  of  homopolymerization  of  length-  and  functionality-distributed 
chains  analyzed  by  Miller  and  Macosko^  using  standard  branching  process  models.  However,  that 
approach  assumes  that  the  characteristics  of  two  chains  with  a  common  crosslink  are  statistically 
independent.  Our  system  of  living  chains  does  not  satisfy  that  assumption.  For  example,  pick  a  random 
crosslink  and  look  at  the  two  chains  joined  at  the  crosslink;  let  their  lengths  be  and  L^.  We  know  that 
the  number  distribution  of  Lj  (and  L2)  is 

L  =  1  +  Ng  +  Ng°  +  N*  +  N“  (56) 


where  Nf ,  Njj°,  N48,  and  N|8  are  independent  Poisson  random  variables  with  means  p|,  pf  °,  p48,  and 
p88  respectively.  Consider  the  event  {  L  =  1  }:  in  order  to  see  a  chain  of  length  1  with  a  crosslink  it  is 

necessary  that  we  have  chosen  a  crosslink  with  an  "i"  unit  and  no  growth  has  occured: 


P(L  =  1)  =  P  (site  =  "i")  P  (Ng  +  Ng°  +  N*  +  Ngg  =  0) 


(57) 


-  2  Kjt  *  Kgit 

2[A“]t  +  2[A‘g]t  +  2[Af]t 
where  p  =  pf  +  p|°+  p4g+  p|g 


exp  ( -p ) 


P(L1  =  1,L2  =  I)  =  P  (both  sites  =  "i" )  P  (no  growth  on  either) 


2[A“],  +  2(A'4‘l,  +  2(A“]1 


[exp  (-H)]2 


(58) 
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For  Lt  and  L2  to  be  independent  it  is  necessary  that 
P(L1  =  1,L2=1)  =P(L1  =  l)p(L2  =1) 

which  is  equivalent  to 


(59) 


Kit 


f 


[A“](  +  [A‘4g]t  +  [Agg]t 


2[A“]t  +  [A‘4g]t 


A 


[A4U],  +  [A‘4g]t  +  [Agg]t 


V 


(60) 


) 


which  is  generally  not  true.  So  the  crosslinked  chains  are  not  statistically  independent  Assuming 
independence  is  probably  not  a  bad  approximation,  especially  when  the  mean  chain  length  is  not  too 
short.  However,  in  this  paper  we  want  to  present  an  exact  model,  subject  to  our  assumption  of  no 
intramolecular  reactions.  Thus  our  crosslinking  model  has  three  types  of  crosslinks. 

Our  crosslink  model  consists  of  chains  with  labelled  chain  segments,  "i"  and  "g".  Crosslinks 
may  be  "ii",  "ig",  or  "gg".  Conditional  on  the  type  of  crosslink,  the  statistical  properties  of  the 
crosslinked  chains  are  independent  This  seems  to  be  an  accurate  description  of  our  situation.  We  can 
analyze  it  combining  ideas  from  our  work  on  crosslinking  chains^  and  our  work  on  nonideal 
systems^  1*12.  We  shall  compute  weight  average  molecular  weight  in  the  pre-gel  region  and  weight 
fraction  soluble  in  the  post-gel  region.  The  number  average  molecular  weight  was  computed  in  eq.  38. 

Weight  Average  Molecular  Weight 

To  compute  the  weight  average  molecular  weight,  we  must  first  compute  the  weight  seen 
looking  from  a  crosslink  ima  a  chain.  There  are  two  possibilities:  we  may  be  looking  into  an  ”i” 
segment  or  into  a  "g”  segment.  We  denote  total  weights  as  E  (W^1)  and  E  (W“),  respectively. 

First  consider  an  "i"  segment.  The  weight  seen  looking  into  this  segment  includes:  mass  of 
root  segment,  mass  of  the  rest  of  the  chain  which  consists  solely  of  "g"  segments,  mass  looking  out  from 
Agg  units  and  A4g  units  on  the  chain.  The  mass  of  the  root  segment  equals  Mj  +  MA4/2.  The  expected 
mass  of  the  rest  of  the  chain  equals  En  (MXGS)  by  previous  reasoning.  The  expected  mass  looking  out 
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from  A88  units  equals  the  expected  number  of  such  units  times  the  expected  weight  looking  into  random 
"g"  segments.  The  expected  mass  looking  out  from  units  is  determined  analogously.  Combining  all 
this  gives 

E(W^ )  =  Mj  +  MA4  /  2  +  En(MXGS)  +  ji“  E(W*  )  +  E  (W”  )  (61) 

Now  consider  a  "g"  segment.  The  weight  seen  looking  in  is  the  same  as  looking  into  an  "i" 

segment,  expect  that  the  root  segment  is  different  and  therefore  the  initiated  chain  segment  must  be  treated 

separately.  We  shall  see  the  mass  of  the  initiating  segment  as  well  as  the  weight  attached  to  it  if  it  is  an 
A^8  or  A^\  This  gives 

E(yT)  =  MAy2  +En(MXGS)  +  jifECWJ1)  +  ^8E(wf) 

+  E  (Mia)  +  P(IA  =  A’8)  E  (W“)  +  P  (IA  =  A“)  E  (W*  )  (62) 


Solving  the  above  two  equations  we  get 

Mj  +  MA4  /  2  +  En  (MXGS)  +  ^8g  E  (Wlgn) 

E  (W*1)  = - - -  (63) 

1  -  H4' 

[1  +  P(IA  =  A“)][M./2  +  E  (MXGS)]  +  (1  - *ij) E(MJ  +  (p^8  +  P(IA=A'i)) M. 

E  (W  m  )= - 4- - : - - - (64) 

(1  -  pj)(  1  -  P(IA  =  A;8))  -  p88(l  +  P(IA  =  A^)) 

The  equation  for  E  (W“)  can  be  computed  and  the  value  plugged  into  the  equation  for  E  (W*1).  Gelation 

occurs  when  the  above  denominator  goes  non  positive,  i.e.  condition  for  gelation: 

(1-pf  )(l-Pa=A‘8))  -  pf  (1  +P(I  =  A“ ))  <S  0  (65) 

In  order  to  compute  the  weight  average  molecular  weight  we  randomly  select  a  unit  of  mass, 
designate  the  monomer  obtained  as  the  root,  R,  and  look  at  the  mass  of  the  root  plus  all  attached  mass. 

=  S  E(W|R)  P(R) 
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=  P(R  =  A^  )M^  +  P(R  =  A12)(MA2  +  EW“-Ma</2) 

+  P(R  =  A*)(MA2  +  EW^'  -  Ma</2) 

+  P  (R  =  A°°)  +  P  (R  =  A;0)  (M^/2  +  EW*n)  (66) 

+  P  (R  =  Af)  (MaJ2  +  EW^  )  +  P  (R  =  A*  )  (EWf  +  EW^  ) 

+  P(R  =  A“)(2EW")  +  P(R=A“  H2EWJ1) 

The  weight  average  of  the  reacted  species  can  be  computed  using  the  probabilities  P’  (R  =  *)  defined 
above. 

Weight  Fraction  Soluble 

To  compute  weight  fraction  soluble  we  pick  a  root  by  mass  and  look  out  from  it  and  ask  if  only 
a  finite  amount  of  mass  is  seen.  Looking  out,  there  are  three  possibilities:  we  see  nothing,  we  look  ime 
an  "i"  segment,  or  we  look  into  a  "g"  segment: 


Define  the  events 

F?n  s  finite  mass  is  seen  looking  into  "i"  segment 
F^n  =  finite  mass  is  seen  looking  into  "g"  segment 
and  let  P  (F*")  and  P  (F^n)  be  the  probabilities  of  these  events  occurring.  Then 
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weight  fraction  soluble  (wfs)  = 

P  (R  =  A°  )  •  1  +  P  (R  =  Aj  )  P  (F*"  )  +  P  (R  =  A*  )  P  (F*gn  ) 

+  P  (R  =  A°°  )  •  1  +  P  (R  =  A^°  )  P  (pi"  )  +P(R  =  A“)P(Fi.n)2  (67) 

+  P(R  =  Af  )P(F^)  +  P(R  =  A*  )P(lf  )P(^)  +  P(R  =  Af  )  P  (F^  )2 

If  one  desires  to  compute  the  weight  fraction  soluble  of  the  yield,  then  he  should  use  the  above  formula 
with  the  probabilities,  P'  (R  =  *). 

Thus  it  suffices  to  compute  P  (Fj®)  and  P  (Fg“).  First  consider  looking  into  an  "i"  segment: 


We  are  interested  in  the  number  of  "g"  segments  on  the  chain  which  belong  to  A*}8  and  Afg  species. 

Recall  that  the  number  of  these  species  are  L^g  and  L|*  respectively,  independent  Poisson  random 

variables  with  means  p^g  and  |ijg,  respectively.  Thus  the  event  Ff1  is  identical  to  the  event  that  all  A* 

species  on  the  chain  lead  to  finite  mass  and  all  A*8  species  on  the  chain  lead  to  finite  mass.  But  the  event 
that  an  A^g  species  leads  to  finite  mass  is  simply  Fjm,  and  the  event  that  an  A8*  species  leads  to  finite  mass 

is  F‘n.  Putting  this  together  we  get 

P(Ff)  -  (£  P  (Ff  )'  P  (I*  =  x))  (£  P  (F™  )*  P  (L«  =  x  )) 

x  =0  X  =  0 

-  exp  (p*  (  1  -  P  (F‘"))jexp  (V4‘  (  1  -  P  (F^) j 
=  *  (P  (F*")  ,  P  (Fi")'| 


(68) 
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Now  consider  P  (F^1).  This  is  more  complicated  because  of  the  initiated  unit: 


The  event  F^1  is  equivalent  to  the  event  that  all  additional  growth  units  on  the  chain  lead  to  finite  mass  and 
the  initiator  unit  of  the  chain  leads  to  finite  mass,  F"ut. 

P(F°ut)  =  P(IA  =  A° )  •  1  +  P(IA  =  aJ0  )  •  1 

+  P(IA-aJ)-  P(Ff)  +  P(IA  =  a|}8)P(F^)  (69) 

Since  growth  segments  of  the  chain  additional  to  the  randomly  chosen  root  segment  of  the  chain  have  the 

same  statistical  properties  regardless  of  whether  the  root  is  an  "i"  or  "g"  segment,  the  event  of  these 
leading  to  finite  mass  is  statistically  the  same  as  F^.  Therefore 

P(F‘gn)  =  P(F;n)  P(F7)  (70) 

Solving  we  get 


P(Fgn )  = 


P(F*™)  J” P(I A  =  A'J)  P(Fi.n)  +  P(IA  =  A°  )  +  P(IA  =  A*4° )' 


1-P(F‘"  )  P(I  =  A'g) 


=  V  (P(F-n )) 


(71) 


The  equation  (68)  becomes 

p(f{")  =  »(p(Fj"),  v(p(F‘n»)  (72) 

and  can  be  solved  numerically  for  the  value  of  P(F^n)  by  the  method  of  binary  search.  The  value  can  then 
be  substituted  into  the  equation 
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P(F^)=v(P(F'i"))  (73) 

to  obtain  the  value  of  P(F“).  These  values  are  then  used  to  compute  wfs  by  equation  67.  This  approach 

could  also  be  extended  to  calculation  of  crosslink  density  following  ref.  3. 

Some  Calculations  and  Comparison  to  Experiment 

The  equations  derived  above  were  solved  using  a  Fortran  program  to  give  Mn  (eq.  38);  Mw 
(eq.  66)  and  wsoi  (eq.  67)  for  the  living  copolymerization  of  A2  with  A4  with  unequal  reactivity.  Since 
the  rate  equations  can  be  integrated  analytically  (eq.  34),  the  program  is  really  just  a  sequence  of 
algebraic  substitutions. 

The  required  inputs  are  the  molecular  weights  of  the  monomers  and  initiator,  MA2,  and 
Mi,  the  ratios  of  their  concentrations  [A4]^[A2]0  and  [IMA^q  and  of  their  rates  of  propagation  r2/ri 
and  ttJx\  (see  eq.  33).  Initiation  is  assumed  to  be  much  faster  than  these  propagation  rates  and 
independent  of  monomer  type. 

Figure  2  shows  typical  results  for  the  simplest  nonideal  case,  homopolymerization  of  divinyl 
monomer,  A4.  The  conditions  were  taken  from  Straehle  and  Funke24  for  the  polymerization  of  ethylene 
glycol  dimethacrylate  (EGDMA  5%  in  toluene)  with  n-butyl  lithium.  Both  vinyl  groups  on  the  EGDMA 
should  have  the  same  reactivity  but  after  one  becomes  incorporated  into  the  chain  we  might  expect  the 
pendent  vinyl  to  be  less  reactive  due  to  topological  hindrance  of  the  chain.  This  has  been  noted  in  free- 
radical  polymerization  of  the  same  monomer.16  Figure  2  shows  that  the  yield  of  polymer  at  gel  is 
increased  from  0.29  to  0.73  when  the  pendent  vinyl  reactivity  is  reduced  by  a  factor  of  12.  The  shapes 
of  the  Mw,  Mn  and  w^j  curves  remain  similar.  Polymer  yield  is  plotted  in  these  figures  because  it  is 
easy  to  measure  experimentally  and  it  is  the  commonly  reported  measure  of  reaction  conversion  (rather 
than  conversion  of  vinyl  groups). 

Straehle  and  Funke  report  a  gel  point  at  about  40%  polymer  which  corresponds  to  about  a 
50%  reduction  in  reactivity.  One  would  expect  a  greater  reduction.  It  may  be  that  their  actual  initiator 
concentration  was  lower  than  reported  due  to  deactivation.  The  value  of  [I]0  given  is  rather  high  and 

would  only  lead  to  a  4000  molecular  weight  polymer  if  there  were  no  branching. 


31 


A  more  complete  set  of  data  on  nonlinear  anionic  homopolymerization  is  given  by  Lutz, 
Beinert  and  Rempp26  for  1,3  diisopropenyl  benzene  10%  in  THF  at  -  30°C  initiated  with  1 
phenylethylpotassium.  At  the  low  initiator  level  used  in  their  experiment  the  gel  point  would  be  at  8.6% 
polymer  if  all  vinyls  were  equally  reactive.  We  see  that  experimentally  the  gel  point  occurs  much  later  in 
the  reaction. 

A  reduction  of  pendent  vinyl  reactivity  by  nearly  500  times  fits  the  experimental  Mw  and  Mn 
values  fairly  well.  The  data  was  taken  by  GPC  which  is  surely  too  low  for  these  branched  molecules  at 
high  molecular  weight  Although  the  delay  in  gelation  can  be  explained  by  reduced  pendent  vinyl 
reactivity  it  could  also  be  due  to  cyclization. 

In  the  general  case  of  A2  +  A4  copolymerization  we  can  have  at  least  three  different 
propagation  rates,  ri  for  the  rate  of  addition  of  the  monovinyl,  A2,  r2  for  either  of  the  divinyl  groups  and 
13  for  the  pendent.  Figure  4  shows  the  effect  of  changing  xjjxx  and  xtJx\  on  the  gel  point.  Initial 
conditions  are  from  Lutz  et  al26  and  are  given  in  the  Figure  caption.  About  10%  of  the  monomer  is 
divinyl  and  the  initiator  concentration  is  typical  (it  would  produce  a  linear  polymer  with  Mw  -  40,000). 

Since  the  ideal  gel  point  is  rather  low,  10.7%,  decreasing  the  relative  reactivity  of  the  divinyl 
or  its  pendent  groups  has  a  much  greater  effect  on  gelation  than  increasing  it  When  the  product 
(r2A‘i),(f3/fi)  <  10-3,  the  gel  point  appears  to  reach  a  limiting  value  ~  0.9. 

We  can  use  Figure  4  to  determine  the  best  parameters  to  fit  the  data  of  Lutz  et  al26  on  the  small 
copolymerization  of  a  methyl  styrene  with  1,3  diisopropenyl  benzene  (1,3  DIB),  see  Table  2.  By 
titration  and  NMR  they  found  that  1,3  DIB  goes  into  the  chain  about  twice  as  fast  as  the  a  methyl 
styrene.  Thus  r2/rj  ~  2.  The  gel  point  was  ~80%  polymer.  From  Figure  4  this  would  indicate  rytj  £ 
0.01.  For  the  calculations  in  Table  2  we  used  13/ri  =  0.008  (or  x?Jx2  =  0.004)  and  get  a  reasonable 
match  to  the  molecular  weight  data. 

It  is  interesting  to  note  that  when  1,3  DIB  was  homopolymerized  (Fig.  3)  a  value  of  xtJx2  = 
0.0022  fit  the  data.  This  is  about  half  the  value  for  copolymerization.  It  is  quite  reasonable  that  the 
pendent  vinyls  on  the  copolymer  will  be  more  reactive  since  they  will  be  much  less  crowded,  spaced 
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apart  by  about  five  a  methyl  styrene  units.  Again  this  discussion  assumes  that  there  are  no 
intramolecular  cycles. 

A  complete  set  of  reactivity  ratios  was  reported  by  Worsfold  for  the  copolymerization  of 
styrene  with  /7-divinyl  benzene.29  He  used  UV  spectroscopy  to  measure  six  of  the  nine  possible 
propagation  rates  and  estimated  the  remaining  three.  He  found,  as  we  assumed  in  our  derivation,  that 

rates  did  not  depend  on  the  nature  of  the  chain  end  but  do  depend  on  the  adding  monomer.  The  rate 
constant  for /7-DVB  adding  to  a  styrene  anion  is  12.3  mol/£«s,  quite  close  to  that  for  adding  to  a  DVB 

anion,  which  is  1 1.2.  However,  the  rate  constant  for  adding  styrene  to  the  chain  end  is  1.16  or  1.12 
mol/ft«s  respectively.  Thus  p-DVB  adds  10  times  faster  to  a  living  chain  than  styrene  monomer.  The 

pendent  vinyl  group  has  nearly  the  same  reactivity  as  the  styrene  0.9  mol/£«s. 

We  can  use  these  data  to  determine  n,  T2  and  rj  and  with  Worsfold's  reactant  concentrations 
we  have  all  the  information  needed  to  calculate  the  gel  point.  From  eq.  65  we  obtain  8.7%  polymer  at 
gel  while  Worsfold  reports  47% !  This  value  was  for  an  initial  monomer  concentration  of  only  8.2%  in 
benzene.  When  he  reduced  the  concentration  to  4.3%  the  gel  point  increased  to  77%  polymer.  This 
effect  of  dilution  on  gel  point  is  a  clear  indication  of  intramolecular  cycles.  Such  large  changes  in 
molecular  weight  and  gel  point  at  high  dilution  are  also  observed  in  stepwise  polymerizations.14'32 

If  we  plot  Worsfold's  gel  points  versus  reciprocal  monomer  concentration  in  the  manner  first 
done  by  Jacobson  and  Stockmayer32  we  get  an  intercept  at  infinite  concentration  in  reasonable  agreement 
with  the  value  we  calculate  (see  Figure  5).  Although  there  is  insufficient  data  to  draw  strong 
conclusions,  this  result  indicates  that  the  rate  constants  are  correct  and  the  main  deviation  is  due  to 
intramolecular  reaction.  Lutz  et  al's  polymerizations  were  run  at  about  10%  monomer  in  THF  so  their 
polymers  must  also  contain  quite  significant  intramolecular  cycles.  The  r3  values  used  to  fit  the  data  in 
Figure  3  and  Table  2  are  probably  much  lower  than  the  true  values,  perhaps  ten  times,  to  compensate  for 
the  additional  effect  of  cyclization. 

Eschwey  and  Burchard19  have  measured  gel  points  for  styrene-divinyl  benzene  systems  at 
2.5%  monomer  in  toluene.  They  found  great  suppression  of  the  gel  point  and  speculate  that  over  90% 
of  the  reacted  pendent  double  bonds  are  involved  in  cyclics. 


33 


Both  Eschwey  and  Burchard  and  Straehle  and  Funke,24  unlike  Lutz  et  al,  report  very  high 
molecular  weights  near  the  gel  point,  a  decade  higher  in  some  cases,  than  would  be  predicted  by  eq.  66 
with  13  adjusted  to  give  their  experimental  gel  point  (see  Figure  2).  This  may  also  be  due  to  extensive 
cyclization  or  perhaps  microgel  and  microheterogeneities  in  the  polymerization. 

CONCLUSIONS 

Following  our  general  paper27  on  treatment  of  combined  kinetic  and  Markovian 
polymerization  problems  we  have  shown  how  the  kinetic  process  of  living,  anionic  polymerization  can 
be  combined  with  the  random  process  of  connecting  these  living  chains  into  a  network.  We  treated  the 
simplest  copolymerization  case:  a  difunctional  plus  a  tetrafunctional,  A2  +  A4  ,  but  the  methodology  can 
readily  be  extended  to  several  monomers  and  other  functionalities.  We  have  modelled  unequal  rates  of 
monomer  addition  but  considered  these  to  be  independent  of  chain  end  type.  Rates  could  be  made 
dependent  on  chain  end  but  this  doesn't  seem  to  be  necessary  for  most  anionic  polymerziations,22*29 
unlike  the  free-radical  case. 

Initiation  rate  was  assumed  fast.  This  is  nearly  always  the  case  in  anionic  polymerization  of 
vinyl  monomers  but  may  not  be  true  for  other  ionic  cases  such  as  epoxy  ring  opening.  In  Table  1  we 

show  that  including  initiation  rate  is  only  critical  for  short  chains.  Variable  initiation  rate  can  be  treated 
by  introducing  another  set  of  differential  equations  for  the  A^,  A°‘.  etc.  species  similar  to  eq.  33. 

Our  calculations  compare  well  to  the  living  homopolymerization  and  copolymerization  data  of 
Lutz,  Beinert  and  Rempp.26  However,  we  had  no  value  for  their  pendent  vinyl  reactivity  and  thus 
adjusted  this  parameter  to  fit  the  data  Reductions  in  pendent  reactivity  of 450  for  the 
homopolymerization  and  250  for  the  copolymerization  were  necessary  for  a  reasonable  fit  to  Lutz  et  al’s 
molecular  weight  results.  Worsfold's  measurement  of  pendent  reactivity  by  UV  spectroscopy  on  a 
similar  monomer  showed  only  a  reduction  of  ten.  We  conclude  that  the  apparent  additional  reduction  is 
reactivity  must  be  due  to  pendent  groups  being  wasted  in  intramolecular  cycles.  Worsfold's  gel  points  at 
different  dilution  support  this  conclusion. 
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Anionic  polymerizations  of  vinyl  monomers  are  typically  run  dilute,  5-20%  monomer,  so 
cyclization  will  be  a  major  factor  in  nonlinear  polymerizations.  Further  experimental  work  like  Figure  5 
is  needed  to  quantify  the  effect  Extension  of  the  theoretical  approach  given  here  to  include  cycles  is  a 
formidable  challenge. 
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Table  1.  Comparison  of  gel  conversion  of  equal  initiation  and  propagation  rates  and  for  very  rapid 
initiation. 


2[A4]o 

[Ho 

ac 

,3 

II 

=  OO 

=  10 

(eq.  16) 

(eq.  30) 

(ref  21) 

2 

.7071 

.8860 

.864 

10 

.3162 

.3317 

.332 

100 

.1000 

.1005 

• 

0 

0 

1000 

.0316 

.0316 

.032 

*  The  actual  value  of  0.105  given  in  ref.  21  was  presumed  to  be  in  error. 


Table  2.  Copolymerization  of  a  methyl  styrene  with  1,3  diisopropenyl  benzene. 

[A4I/IA2] 

in  polymer  Mn  Mw 


yield 

exp. 

ca.  40 

GPC 

eq.  36 

GPC 

eq-  66 

.057 

.36 

.35 

5300 

4200 

5750 

4300 

.31 

.26 

.25 

20,600 

24,200 

24,500 

25,800 

.62 

.17 

.16 

47,900 

62,000 

91,500 

95,700 

from  Lutz  et  a l26  MA2  *  1 18,  MA4  =  158,  Mi  106,  [A^o^fo  =  0.0947,  [A^q/Mo  =  556  (the  effective 
initiator  concentration  ws  reduced  to  2/3  the  reported  value  to  match  the  low  molecular  weight  data); 
calculations  assume  r2/ri  =  2,  ttJt\  =  .008.  Reactant  concentration  was  14%  in  THF. 
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Figure  Captions 

Figure  1.  Molecular  weight  build  up  during  ideal  polymerization  of  a  tetrafunctional  monomer  by  three 
different  mechanisms:  chain  ( - eq.  21),  living  ( - ,  eq.  17)  and  step  ( - ,  eq.  19). 

The  "typical"  conditions  are  :  kinetic  chain  length  of  50  for  the  chain  addition  (q  =  0.98)  and 
2[A4]0/[I]0  =  100  for  the  living  system.  For  the  dashed  lines  to  the  right  these  parameters 

have  been  reduced  to  match  the  gel  point  of  the  stepwise  polymerization  which  will  always 
gel  at  ac  =  l/(f  -  1)  or  1/3. 

Figure  2.  Number  average  and  weight  average  molecular  weight  and  soluble  fraction  vs.  polymer  yield 
(fraction  of  monomer  incorporated  into  polymer).  A4  homopolymerization: - ideal; 

- pendent  vinyls  are  12  times  less  reactive.  MA4  =  198,  Mj  =  56;  2[A4]q/[I]q  =  43.4,  the 

conditions  reported  by  Straehle  and  Funke  to  polymerize  ethylene  glycol  dimethacTylate. 

Figure  3.  Homopolymerization  of  1,3  diisopropenyl  benzene.  ■  Mw  and  ▲  Mn  data  et  al.26  Open 

points  are  believed  to  be  in  error  (low)  due  to  the  method  used,  GPC.  Solid  curves  are  Mw 
from  eq.  66  at  various  ratios  of  pendent  vinyl  to  monomer  reactivity  ryr2  .  Dashed  curve  is 
for  Mn  (eq.  36)  for  r3  /r2  =  0.0022. 

Figure  4.  The  effect  of  unequal  reactivity  on  yield  at  gel.  The  ideal  gel  point  occurs  at  0.107  polymer 

fraction.  Decreasing  the  reactivity  of  the  divinyl  monomer,  r2,  or  the  vinyl  pendent  on  the 
chain,  r3  ,  with  respect  to  the  monovinyl  monomer,  rj ,  increases  dramatically  the  fraction  of 
polymer  at  gel.  [A4]o/[A2]0  =  0.0974;  [A2]o/[I]0  =  556. 
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Figure  5.  Gel  point  vs.  reciprocal  monomer  concentration  for  copolymerization  of  styrene  with  p- 
divinyl  benzene. 

[A2^  =  0.374  moW,  [A^/ [A2]Q=  4.28  x  10*3,  ^/IUq  =  4560 

□  ^  =  0.714 moW, [A4]0/[A2]0=  1.57  x  lO'3^  /[T]Q  =  8500 

O  calculated  for  same  conditions  as  0  with  r,  =  1.14,  r2  =117, 

29 

r3  =  0.9,  the  average  values  given  by  Worsfold. 
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